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Abstract. 

We consider a population of N labeled random walkers moving on a substrate, 
and an excitation jumping among the walkers upon contact. The label X(t) of the 
walker carrying the excitation at time t can be viewed as a stochastic process, where 
the transition probabilities are a stochastic process themselves. Upon mapping onto 
two simpler processes, the quantities characterizing X(t) can be calculated in the 
limit of long times and low walkers density. The results are compared with numerical 
simulations. Several different topologies for the substrate underlying diffusion are 
considered. 
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1. Introduction 

A general stochastic process can be viewed as the time evolution of one (or more) 
random variable pQ, the particular dependence on t of the transition probabilities 
between the states giving rise to different models. Among the most widely studied 
stochastic processes in physics are Markov processes, where the transition probabilities 
at ti > t depend only on £(£) and t, and not on the previous history of the system. In 
the simplest case the time parameter t is discrete, and is called a Markov chain; 
the case of a Markov chain with transition probabilities independent of t is by far 
the most studied. If the transition probabilities in the time interval (to, to + t) do 
depend on t (with a given distribution function), but not on t 0) we have homogeneous 
processes. Depending on the particular functional dependence on t, we can obtain 
Poisson processes, Wiener processes, and so on. Relaxing the homogeneity property, we 
can obtain the inhomogeneous version of the previous processes. 

Much more general assumptions on the time- dependence of the transition 
probabilities can be given, but the resulting models are rarely explicitly solvable. In this 
paper we define and solve a particular discrete-time stochastic process: its transition 
probabilities are a stochastic process themselves. 

The process we consider is a "second-level" random walk, or random walk on 
random walkers. We consider N labeled random walkers, diffusing on a given substrate. 
Such random walkers can define a dynamic meta-graph: each random walk is seen as 
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a node of the meta-graph and a link between two of them is drawn whenever they are 
within a distance R on the substrate. Then, we study the diffusion of a "second-level" 
random walk on such meta-graph. 

Apart from its mathematical interest, this kind of system is also able to model a 
diffusion-reaction process. In fact, each walker diffusing on the substrate represents 
a particle (all particles belonging to the same chemical species) that can be either in an 
excited (A*) or in an unexcited (A) state, the former corresponding to the node carrying 
the second-level random walker. When an excited particle meets an unexcited one, they 
immediately react according to the scheme 

A* + A^A + A*. (1) 

This reaction mechanism is known as homogeneous energy transfer (ET) which takes 
place from an excited molecule [donor (A*)] to another unexcited molecule [acceptor 
(A)], according to the scheme (CQ). This process stems from Coulombic (long-range 
[2]) and exchange (non-radiative, short-range [Hj) interactions amongst the particles. If 
we just focus on the energy transfer via exchange (under the implicit assumption that 
the relaxation takes zero time), this allows to restrict transfer interaction to nearest- 
neighbour particles only. 

If we define an abstract space whose points are the N random walkers, the excitation 
transfer corresponds to a stochastic process X(t) on the points of this space; hence, 
to a "second-level" random walk. The transition probabilities of this process depend 
on the relative positions of the random walkers, hence they are a stochastic process 
themselves. It is possible to show how the process X(t) can be mapped exactly onto 
simpler processes, involving N or iV — 1 simple random walkers on the same lattice; the 
study of the excitation jumps is here mapped on the study of the passage times of these 
walkers through the origin. These simpler processes can be solved in the limit of large 
times and low walkers densities. 

The paper is organized as follows. In Sec. [2] we describe the model; in Sec. [3] 
we provide two mappings to simpler processes that allow us to obtain the asymptotic 
behaviour of the quantities of X(t). In Sec. H] these results are compared with numerical 
simulations. Sec. [5] contains our conclusions and perspectives. 

2. The model 

We consider N regular random walkers, labeled with the numbers from 1 to N, moving 
on a finite structure (henceforth, the substrate). The position of the ith walker at time 
t is Xi(t); at time all the positions are random. At t — one of the walkers, io, carries 
an excitation; we assume without loss of generality that Zq — 1. 

The following usual quantities for random walks on lattices will be useful. For a 
walker starting from r at time 0, we define the probability Po(r, t) of being at at time 
t, and the probability Fq(v, t) of being at for the first time at time t. We also define 
their generating functions, Po(r, A) = X^o-^o( r ; fyX* and F (r, A) = Y^i^o( r i 
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We fix a collision radius R > 0: at time t two walkers meet (or collide) if their 
distance on the lattice is < R. In this paper we consider R = 1 but there are no 
substantial differences for different R (the choice R = is here neglected to avoid parity 
effects, and used for explanations only in Sec. [3]). When the walker % carrying the 
excitation collides with another walker j, the excitation jumps from i to j. If it collides 
with more than one walker at the same time (which we will call a multiple hit), the 
excitation jumps on one of them chosen randomly. 

The model just described defines a discrete-time stochastic process X(t), where the 
state space of the system is composed by the set of the random walkers. At time t the 
system is in state % if the excitation is on walker i 

Formally, the process is defined by the state space 

• X(t) e AT, AT = {1,2,..., N}; 

by the initial condition: 

• X(0) = 1; 

and the evolution rule: 

• let X(t) = i; consider the set C = {j : \\xi(t) — Xj(t)\\ < R; j ^ i} [j]. 
If C = 0, then X(t + 1) = %. 

If C 0, then X(t + 1) = j, where j is chosen randomly among the elements of C 
with equal probability. 

Here, the transition (or jump) probabilities, given by the evolution rule, are a 
stochastic process. In particular, at time t the transition probability from state i 
(X(t) = i) to state j (X(t + 1) = j) is a function of the positions Xi(t) and Xj(t) 
of the two RWs, hence a function of two stochastic processes. 

Several quantities can be defined for X(t), much in the same way as for regular 
random walks on a lattice. We define: 

• J(t), the average number of jumps performed by the system up to time £; the 
probability J~(t, h) that the number of jumps performed by the system is h at time 
t, J(t) = Eti kj(t,k). 

• S(t), the average number of different states visited at time t; the probability 
S(t,k) that k different states have been visited by the system at time t, S(t) = 

• the Cover Time r, defined as the average time required to visit all the N walkers 
(analogous to the lattice-covering time for random walks [7]). We also define n as 
the average number of jumps required to visit all the states (n < r). 

The substrates considered will be Euclidean (hypercubic) lattices of linear size L 
and volume L d (with d — 1, 2, 3), endowed with periodic boundary conditions. 

X Here, \\x — y\\ denotes the chemical distance between x and y, both for Euclidean and fractal lattices. 
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We also will consider fractal substrates. It is well known [U |5] that fractals are 
described by at least two different dimensional parameters. One is the fractal dimension 
dp, describing the large-scale dependence of the volume (or mass) V(r) of the structure 
on the distance r from a point chosen as the origin: Vo(r) ~ aor dp (here and in 
the following lines, a , b and c are constants depending on the point 0). The other 
is the spectral, or connectivity, dimension d s , describing the long-time behaviour of 
diffusive phenomena on the fractal. For example, for t — > oo the probability of return 
to the starting point for a RW on the fractal is Pa(t) ~ bot~ ds ^ 2 , and the average 
number of different sites visited by the RW is S(t) ~ cot min< - ds//2 ' 1 \ For Euclidean 
lattices, d s = dp = d. In a lattice (either Euclidean or fractal) with d s < 2 a random 
walker starting from a point is bound to return to an infinite number of times with 
probability 1, and the lattice is called recurrent. For d s > 2, the walker has a non- 
null probability to escape to infinity without returning to 0, and the lattice is called 
transient. 

The fractal lattices we will consider (fig. [1]) are Sierpinski gaskets of linear size L 
and volume L log3//log2 (dp = log 3/ log 2). Their spectral dimension is d s = 2 log 3/ log 5 
(hence, they are recurrent: d s < 2) . 

All the quantities we are interested in will be examined as functions of N and L. 

3. Analytical Results 

The purpose of this section is to show how our model can be mapped onto two different, 
and easier, models, that we shall call picture 1 and 2 respectively. In these two pictures, 
and in the low-density (LD) limit (when multiple hits are negligible), the asymptotic 
behaviour of the quantities of the previous section can be found. 
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Let us take figure [2] as a reference. The upper part of the figure exemplifies the 
basic process. At t — the excitation is on walker 1 (the system in state 1); at t\ walker 
2 hits walker 1 and the excitation jumps on walker 2 (the system jumps on state 2). At 
times t 2 and £3 the excitation jumps on walker 3 and then on walker 1 again. This can 
be summarized by introducing the sequence of jumping times 

0, h, t 2 , t 3 , . . . (2) 

and the sequence of visited states 

0, h, i 2 , 13, ■■■ (3) 

Picture 1 (stuck- and- free picture) We consider the process in the reference frame 
of the excitation. In this frame, the walker carrying the excitation is stuck at the origin, 
and the other N — 1 walkers perform a regular random walk, with 2 jumps on each 
time step. Here, the jump of the excitation from walker i to walker j corresponds to 
the following: walker j hits the origin and gets stuck, while walker i gets free and starts 
performing its own RW. 

In this picture, the process is a double-state RW process [8], because each walker 
can exist in two different states: either stuck at the origin or free. When a walker is 
free, this picture allows us to use well-known quantities from random-walk theory: for 
example, the probability for walker i, starting from at time 0, of getting stuck at 
the origin at time t is (neglecting multiple hits) F Q (ri — r Q , 2t). This problem is still 
completely described by the above sequences of times ([2]) and states ([3]). 

We remark that this mapping is possible only for translationally invariant 
(Euclidean) lattices, where the lattice in the reference frame of the excitation is the 
same as the original one. It is not possible for fractal lattices; this will be clarified 
below. 

Picture 2 (label permutation picture) 

When walker 2 hits walker 1 at the origin and gets stuck (picture 1), the random 
walk subsequently performed by 1 is just the random walk that would have been 
performed by 2 if no sticking effect had existed: that is, if walkers 1 and 2 simply 
had switched their labels without changing their state. This label switch can be seen as 
the action of a transposition (12) of the numbers 1 and 2 on the sequence N '. 

Consider the process (let us call it the associated free process) with N — 1 free RWs, 
labeled from 2 to N, on the same lattice, and walker 1 stuck once and for all at the 
origin. The process in picture 1 is the same as the associated free process, plus the 
following condition: when a walker hits the origin it switches its label with the last 
walker that has hit the origin before it (with the condition that the first walker has been 
1). In general, when walker i of the associated free process hits the origin, a permutation 
IT = (1 i) of elements 1 and i is induced on the original sequence M (since the last stuck 
walker is always at the first place in the permutated sequence). 

The sequence of jump times ([2]) hence is equal to the sequence of crossing times of 
N — 1 random walkers through the origin. Hence, the sequence of the walkers that cross 
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Figure 2. Top: the original process on a square lattice at four nonconsecutive times 
0, t\, t2, t%. The walker carrying the excitation is the black circle. The excitation 
jumps from 1 to 2 at time t\, from 2 to 3 at time and from 3 to 1 at time £3. 
Middle: picture 1. The same process in the reference frame of the excitation (small 
black circle fixed at the origin). The walkers are stuck at the origin when carrying the 
excitation in the original model, and get free when the excitation jumps to another 
walker. The jumping times are the same. Bottom: picture 2. Here, the black circle 
marks the origin. The associated free process, with N — 1 random walkers labeled with 
the numbers from 2 to N, is shown. Picture 1 is obtained as follows. We start from 
the ordering (1, 2, 3). Each time walker i of the associated process crosses the origin, 
it exchanges its label with the previous walker that crossed it, starting from walker 
1; alternatively spelled, walkers at position 1 and i of the present ordering exchange 
their labels. Hence, 2 crosses the origin at time t\ and exchanges its label with 1; the 
new ordering is (2, 1, 3); 3 crosses the origin at time t 2 and exchanges the label with 
2; the new ordering is (3, 1, 2). Finally, 2 crosses the origin at time t 3 ; the walkers at 
positions 2 and 1 exchange their labels: the new ordering is (1, 3, 2). 



the origin in the associated free process 

0, jl, fa, J 3 , 

is related to the sequence for the original process by 

h = (ni./^; i2 = (n 2 A0 i2 ; is = (n 3 A0j 3 ; • • • 

where III = (1 ji); U 2 = (1 ji)(l j 2 )\ n 3 = (1 ji)(l j 2 ){l j 3 ), and so on. 

Two observations are necessary at this point. First: both pictures are valid only 
for translationally invariant lattices; for fractals, for example, the lattice in the frame 
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of reference of the excitation does not coincide with the original one (indeed, it is not 
even fixed but changes with t). However, several of our numerical results suggest that 
the asymptotic results derived in the euclidean case also hold (in some averaged sense) 
for non-integer-dimensional cases. This point will be stressed again case by case. 

Second: we depicted pictures 1 and 2 for a model with null range R = 0, while most 
of our numerical result concern the case R ^ (mostly R — 1), chosen to avoid the 
parity effects (since most of our lattices are bipartite graphs, walkers starting from the 
"wrong" sites would never meet). A non-null range in the original model corresponds 
to a sticking area greater than the origin in picture 1, and to the passage to a region 
greater than the origin in picture 2. This means that in picture 1 and 2 the walkers can 
perform jumps to the origin even when the origin is not a nearest-neighbor site. We 
expect, however, that the existence of a non-null range will only result in a rescaling 
of the asymptotic laws (usually by a factor v/V, where v is the discrete volume of the 
region). We will stress this point in the analytic results where necessary. 



3.1. Number of jumps for large times 

This quantity is easily calculated in picture 2. If we consider low-density systems, that 
is, we neglect the probability of multiple hits of the origin by the walkers, the number 
of jumps at time t is the number of passages through the origin made by iV — 1 RWs at 
time t, that is N — 1 times the number of passages through the origin made by a single 
RW. The mean number of times that a RW starting from r visits the origin in a walk 
of t steps is independent of r for large t, and equals ~ y, where V is the volume of the 
lattice [9]. The average number of jumps is given by the mean number of times that 
iV — 1 independent RWs hit the origin, that is 

N - 1 

J(t) ~ -^t, (4) 

neglecting multiple hits. In the case of walkers with non-null radius of action we must 
consider a finite-size trap. If v is the volume of the trap, the result is 

J(t) ~ v t- (5) 

For example, for a radius R = 1 we have v = 2d + 1 for hypercubic lattices of dimension 
d. 

For J~(h, t) (the probability that the number of passages performed by the excitation 
is h at time t) no analytical results are known, and we will rely only on numerical 
simulations. 



3.2. Cover Time 

The Cover Time is defined as the average time needed for the system to visit all the 
states. In the LD limit this is equal (looking at picture 2) to the time needed for N — 1 
different walkers to be absorbed into a trap located at the origin. This is a many-body 
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problem (already formulated in the frame of extreme value statistics, see e.g. [B]), and 
its exact solution is not yet known. 

We will adopt here an approximation. We recall that Fq(v, t) is the probability 
density for the first-passage time to the origin of a walker starting from r. We know 
that on hypercubic lattices the average first passage time for a RW through the origin, 
averaged over all possible starting positions, is 



1 °° 

® v ^EvE^'*) ~ a d g d (V), 



r t=Q 

where the approximation is valid for V large; a d is a constant that depends only on d, 
and g d (V) is the volume-depending part: 

V 2 d=l 

9d(V)= { VlogV d = 2 (6) 
V d>2 

In the case of fractal lattices, the general formula (t)v ~ a ds 9d s {V) can be heuristically 
justified, and has been calculated analitically in two particular cases [TUHH]; here, 

V 2 ' d ° d s <2 

9d 8 {V) = { VlogV 4 = 2 (7) 

V d s >2, 

d s being the spectral dimension of the lattice. 

Our approximation consists in assuming that the first passage time of the first 
out of m RWs is that of one RW divided by m. Hence, the time of absorption of the 
first walker is g d (V)/(N — 1), that of the second walker (the first out of iV — 2 left) is 
g d (V)/(N — 2) and so on. The Cover Time is: 

N-l 

t(N,V) ~ ~ b + lo § N + OiN- 1 )} a d9d (V), (8) 

n=l 

where the last relation holds in the limit of large N. 

From what said before, we can easily estimate the average number of jumps required 
to visit all the states: 

N-l 

7r (N,V) = —^r(N,V). (9) 

In fact, as stated by equation (jlj), the average time taken by the excited particle to meet 
another particle out of the remaining iV — 1 is just jfzj- 

3.3. S(t), number of distinct particles visited at time t 

In the low-density limit (again looking at picture 2), this quantity is the average number 
of particles (out of N — 1) that survive at time t with a trap in the origin. This in turn 
is iV — 1 times the survival probability of a single walker with a trap in the origin. 
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This quantity has been calculated in [T2] for Euclidean lattices; let us quote here 
the main results. Let U(t) and S(t) be the survival probability of the walker and the 
average number of sites visited by the walker at time t, respectively. The two quantities 
are related by the formula U(t) — 1 — S(t)/V. Let 5(A) be the generating function of 
S(t) with respect to time. We have S(X) = /(A)/(l — A), where 

/(A)= [(1-A)0(O,A) + 1/L d ]- 1 . 

The function 0(0, A) constitutes the non-singular contribution to the generating function 
P (0, A) as A — > 1. More precisely, 0(0, A) is just a finite sum of terms involving the 
structure function of the substrate. 

The behavior of /(A) near its radius of convergence is governed by 0(0, A), where 
A is the root with the smallest magnitude of the equation /(A) -1 = 0. For d — 1 this 
value is known exactly to be 0(0, A) = 2L/tt 2 . For d — 2, it is found numerically that 
0(0, A) ~ 0.441ogL. For d = 3, A = 1 and 0(0, A) = 1.51.... Given these results, the 
behavior of U(t) for large times is 

U(t) ~ exp f- * ) . (10) 



^0(0, A) 

We will find it expedient to write U(t) ~ e - x dt/g d (V) ^ c f r _ e q Ua ti ns (EI) and (JTj)), 
where all the constants are absorbed in A^. Hence, 



S(t) ~ (N 



cxp 



(11) 



9d(V)J\ ' 

Now, by comparing S(t) ~ V [1 — U(t)] with S(t) we can derive that the fraction of 
distinct particles excited just corresponds to the fraction of distinct sites visited 
by a regular random walker on the substrate. Equation fllip holds also for fractals, 
replacing d with d s . 

For earlier times, the role of topology in the behavior emerges [T2] : 



\ j.inin(d a / 2,1) ' 

U(t) ~ exp ( — — 7T7\ ) • (12) 



9d s {V) 

Finally, notice that the (finite) size R of the trap does not qualitatively affect the 
previous relations while, in general, the value of the constant A^ s may non-trivially 
depend on R. We will deepen this point later in Section 14.21 



3. 4- S(k,t), probability distribution function for the k distinct agents visited at time t 

S(k, t) corresponds, in picture 2, to the probability that the number of walkers absorbed 
into a trap at the origin is k. Recalling that U(t) is the probability that a given walker 
has survived up to t, we have: 

s(k,t) = u(t) N -\i-u(t)) k - 1 1 v ; i. i<k<\ 




Random walk on a population of random walkers 



10 



that is (recalling that for Euclidean lattices d s = d): 

S(k,t) = e -(N-i)x ds t/ 9ds (v) ( e x ds t/ 9ds (v) _ ^fc-i f N-l \ (13) 

Notice that, in the thermodynamic limit, equation ffl3|) becomes a Poissonian 
distribution with average /i = \d s (N — l)t/gd s (V) (see Fig. Cj). 

The time t P eak(&), each distribution is peaked at, can be directly derived from 
equation (fl"3l) : 

V f N - 1\ 

t peak (k) = -log — — . (14) 



A d °\N-k, 

An important feature concerning S(k, t pea k(k)) is that it exhibits a minimum for 
k = k = as can be deduced from equations f|T3|) and (1X41) . 

It is as well possible to calculate the average time r^-k spent by the system having 
visited exactly k different states: 

oo 

TK _^g s(MH __, (15) 

where the last relation was derived in the continuum limit for t. 



4. Numerical Results 



We first consider final quantities, i.e. quantities measured when the excitation has 
covered the whole population of walkers. Subsequently, we will take into account the 
temporal evolution of the system by discussing quantities such as the average number 
of distinct walkers S(t) visited at least once by the excitation, as well as S(t,k) and 
J{t, k) representing the probability distribution of having k distinct walkers visited at 
time t and of having h jumps performed at time t. 




10 1 10 z 10 3 10 1 10 z 10 3 

N N 



Figure 3. Rescaled Cover Time t(N, V) versus the number of walkers making up the 
system and diffusing on a periodic chain (left panel) and cubic lattice (right panel). 
Different sizes are considered, as shown by the legend. Equation ([8]) provides the best 
fit when reactants concentration is small. 
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4-1. Cover Time and Cover Jumps 

In this section we focus on numerical results concerning the Cover Time r and the Cover 
Jumps 7r. We recall that r has been defined as the average time it takes the excitation to 
reach all the N walkers diffusing on the substrate considered. Analogously, ir represents 
the average number of jumps performed by the excitation within the time at which 
S = N. Obviously, n < r. 

In Figs. [3] and H] a proper rescaling of data points confirms the analytical results 
discussed in the previous section (see equations (|Sj) and ©). In particular, in the low- 
density regime, t(N, L) and tt(N, L) depend separately on N and L and their functional 
form is strongly affected by the topology of the lattice underlying the propagation (for 
example notice that for transient substrates 7r gets independent of the size of the lattice). 

4-2. Distinct walkers Visited 

In Section Q] we introduced S it) as the average number of distinct walkers which have 
been excited at least once at time t. 

In Sec. [3] we analytically showed that, in the long-time regime, independently of 
the (finite) substrate topology S(t) grows exponentially with time (see equation ffTTj) ). 
On the other hand, in the early-time regime, for recurrent substrates, a functional 
dependence on the topology is expected, consistently with what found for a random 
walker on a finite lattice [12] . 

Let us first consider the case of a cubic structure for which the behavior of S(t) 
is not expected to display any crossover in time. Indeed, Fig. [5] confirms this: on the 
whole range of time, equation ( TTTj) is a good estimate for S(t) when the density is low. 




Figure 4. Rescaled Cover Jumps n(N, V) for a system of walkers diffusing and reacting 
on a Sierpinski gasket (left panel) and on a cubic lattice with periodic boundary 
conditions (right panel). Different sizes are depicted, as shown by the legend. Equation 
<j9j) provides the best fit when the reactants concentration is small. Notice that in the 
latter case n(N, V) is independent of V. 
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The slope of Vlog ^1 — also allows to derive an estimate for the constant A^ s . By 
fitting numerical data we find that A 3 ~ 2.65(5), A 2 ^ 6.84(4), Ai ~ 10.01(8), (to be 
compared with those in Sec. 13.31 recalling that here R — 1). 

Now, let us consider low- dimensional substrates. The numerical simulations 




Figure 5. Rescaled number of distinct particles visited by the second-level random 
walker as a function of time for a periodic cubic substrate. Equation (llip holds for 
any (low) concentration chosen. The only free parameter in the fitting procedure is 
A^* = 2.65 ±0.05. 




Figure 6. Time dependence for the number of distinct walkers excited at least once 
and diffusing on Sierpinski gaskets of different generations, as shown by the legend. 
The crossover between the two time regimes is apparent by comparing the plot in the 
large figure and the set of data depicted in the inset. The quantity — log(l — S(t) /N)V 
scales respectively as t and t ds / 2 . 



Random walk on a population of random walkers 



13 



0.8 



0.9 




0.7 



-3 












2000 



4000 



6000 



8000 



10000 



12000 



Figure 7. Main Figure: Probability distribution S(k,t) versus time t for a system 
of N — 32 walkers diffusing on a cubic lattice sized L = 16 with periodic boundary 
conditions. Each curve represents a different (even) value of k: starting from the 
leftmost distribution k — 2, 4, 6, 32. Inset: Probability distribution S(k,t) versus 
number of visited random walkers k; three different instant of time are depicted in 
different colors: t — 6 • 10 2 , 1.1 • 10 3 , 2 • 10 3 . Data points (o) are fitted by a Poissonian 
distribution with average [id s — ^d s pt m agreement with what stated in Sec. 13.41 



performed on the chain and on the Sierpinski gasket (see Fig. [6]) support what previously 
stated. In particular, for the latter we show that, at long time, S(t) increases 
exponentially, analogously to what previously found for the cubic lattice. Conversely, 
at small times, deviations emerge: the pure-exponential growth is replaced by e tdsl2 in 
agreement with equation (TT2|) . 

In Sec. [I] we introduced the function S(k,t), representing the probability that, at 
time t, the number of walkers visited at least once by the excitation is S(t) = k. In Sec. [3] 
we also derived a mean-field approximation for this quantity, valid in the low-density 
regime. We now discuss the pertaining results from numerical simulations. 

In Fig. [7| the probability distribution S(k,t) is fitted by a Poissonian law with 
average \x linearly dependent on the density p = y of the system. Moreover, the time 
tpeak each distribution is peaked at depends on k and it diverges logarithmically when 
k — > N (see Fig. according to equation (JHj). 

From the distribution S(k,t) it is also possible to measure the average lifetime 
for the k-th state. This quantity diverges linearly as k — > iV as shown in Fig. [8] where 
results for the cubic lattice are depicted and fitted consistently with equation f fl~5l) . 

An important feature emerging from Fig. [7] is the existence of a minimum for 
<S(k, tp ea k). Indeed, there exists a value k at which the distribution is maximally spread; 
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Figure 8. i pea k(fc) and r N as a function of k for a periodic cubic lattice. The 
dashed lines (whose equations are reported) represent the best fits in agreement with 
equations lfT4|) and (fT5|) . The only free parameter is Xd s and we get X^ 1 = 2.79 ±0.07. 
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Figure 9. Probability distribution J(h, t) versus time t for a system of N — 32 walkers 
diffusing on a cubic lattice sized L = 16, with periodic boundary conditions. Several 
curves are depicted, each referring to a different number of passages h (selected one 
every 5 entries). As h increases, the extremal point of the related distribution t pea k(M 
gets larger, distributions are more and more overlapped and fluctuations get more 
important. The best fit for J(h, t pca k(/i)) is represented by the black line y = A t B , 
with A = 1.11 ± 0.02, B = 0.53 ± 0.01. Data have been averaged over 1.8 • 10 5 
realizations. 



in the average k = y and, correspondently, the statistical knowledge we have about the 
system is minimum. From equation (jTTT) we can estimate t fa ^— log 2. 

Finally, in Fig. [9] numerical results for J{h,t) are depicted. We recall that J{h,t) 
just represents the probability that the number of passages performed by the excitation 
is h at time t. From the perspective of the energy-transfer mechanism this quantity 
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is also of practical interest, especially in the case we allow for energy dissipation or 
emission during transfer. As shown in Fig. [91 there is no extremal point for the envelop 
of such distributions which is indeed characteristic of S(k,t). 

5. Conclusions and perspectives 

We have introduced and studied the diffusion of an excitation (or second-level random 
walker) on a population of N random walkers diffusing on a given lattice (substrate) with 
finite volume V. This results in a stochastic process X(t) whose transition probabilities 
are themselves stochastic. The interest in this kind of problem is also motivated by the 
fact that it provides a model for systems of particles interacting by means of exchange 
energy transfer. 

We showed that in the low-density regime (p = y <C 1) X(t) can be mapped onto 
simpler processes, which allows the analytic calculation of the quantities characterizing 
the diffusion of the second-level RW. This analytic approach becomes rigorous only 
for homogeneous substrates, but yields reliable results also for fractal substrates. We 
presented numerical results supporting our analytical findings. 

There are two main possible developments for this model. First, one can introduce 
a number N e > 1 of excitations jumping among the walkers. This would allow for 
the existence of several donors (excited walkers) in the system at the same time, and, 
possibly, of several excitations residing on the same walker. The rules governing the 
interaction between two donors (i.e., the existence of constraints on the number of 
excitations on a single walker) would have to be included in the model. 

The second development consists in adding more levels of diffusion. If we define a 
set of iV e > 1 excitations, we obtain a set of N e > 1 second-level stochastic processes. 
We can then define a collision rule for those stochastic processes (for example, two of 
them collide when the two excitations are on the same walker). Then, we can introduce 
a third-level stochastic process by allowing a third population of walkers diffuse on the 
second population (that of the excitations). The interplay between the properties of the 
second- and third-level stochastic processes (and a fourth-level one, and so on) could 
then be studied. 
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